Surface and smectic layering transitions in binary mixtures of parallel hard rods 
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The surface phase behavior of binary mixtures of colloidal hard rods in contact with a solid 
substrate (hard wall) is studied, with special emphasis on the region of the phase diagram that 
includes the smectic A phase. The colloidal rods are modelled as hard cylinders of the same diameter 
and different lengths, in the approximation of perfect alignment. A fundamental-measure density 
functional is used to obtain equilibrium density profiles and thermodynamic properties such as 
surface tensions and adsorption coefficients. The bulk phase diagram exhibits nematic-smectic and 
smectic-smectic demixing, with smectic phases having different compositions; in some cases they are 
microfractionated. The calculated surface phase diagram of the wall-nematic interface shows a very 
rich phase behavior, including layering transitions and complete wetting at high pressures, whereby 
an infinitely thick smectic film grows at the wall via an infinite sequence of stepwise first-order 
layering transitions. For lower pressures complete wetting also obtains, but here the smectic film 
grows in a continuous fashion. Finally, at very low pressures, the wall-nematic interface exhibits 
critical adsorption by the smectic phase, due to the second-order character of the bulk nematic- 
smectic transition. 

PACS numbers: 64.70.M-,61.30.Hn,61.20.Gy 



I. INTRODUCTION 

The wetting behavior of molecular smectic (S) liquid 
crystals in contact with a solid substrate [H-Q or at their 
vapor-liquid crystal interfaces (6l4l0| has been an active 
research area since the 80's. Partial or complete wet- 
ting behaviors of the S phase have been found when the 
isotropic (I) or nematic (N) phases are stable at bulk. 
Some liquid crystals exhibit a sequence of first-order lay- 
ering transitions on decreasing the temperature slig htly 
above the IS or NS bulk transition temperatures |l|-[lOf . 
indicating that the partial and complete wetting regimes 
can be mediated by a finite (or infinite) sequence of step- 
wise layer-adsorption transitions. 

Usually the wetting behavior depends on specific in- 
teractions between the surface and the liquid crystal 
molecules, and molecular characteristics such as dipoles 
or length of alkyl chains. Depending on the strength of 
the interactions, solid substrates inducing orientational 
ordering of molecules may favor partial or complete wet- 
ting of the substrate by the N or S phases when the I 
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phase is stable at bulk, while those substrates promoting 
orientational disorder favor partial wetting behavior. 

Freely-suspended smectic films, consisting of a few 
smectic layers surrounded by vapor (V), are formed by 
some liquid crystals and constitute another example of 
phase transitions induced by the presence of a surface. 
These films exhibit the so-called thinning transitions 
whereby the film thickness decreases stepwise as one or 
several layers (depending on the film heating rate) melt 

[IMl. 

As usual in statistical mechanics, lattice models were 
the first to be applied to the study of the surface phase 
behavior of smectic liquid crystals. For example, a ver- 
sion of the Lebwohl-Lasher model, extended to include 
the smectic phase, was used to study the systematics of 
layering phenomena (l3 |. 

Density-functional theory (DFT) has been also suc- 
cessfully applied to the study of the surface phase be- 
havior in liquid crystals adsorbed on solid substrates. 
The extension of the MacMillan theory to non-uniform 
phases with the inclusion of surface interaction potentials 
accounted for layering and thinning transitions (l5l . [l6| . 
However, DFT models that incorporate repulsive inter- 
actions (reflecting molecular volume and shape) using ei- 
ther the Local- (LDA) [l3| of Weighted-Density Approx- 
imation (WDA) [l^, plus anisotropic attractive interac- 
tions via a mean-field approximation, turned out to be 
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more realistic models for the calculation of surface phase 
diagrams. This is due to the fact that (i) the liquid crys- 
tal bulk phase behavior (e.g. values of coexistence den- 
sities and orientational order parameters) is better cal- 
culated from DFT, and (ii) interfacial properties, such 
as the width of the interface or the oscillatory behavior 
of the density profiles, are much better accounted for, 
due to the proper inclusion of pair correlations between 
particles. 

For example, the wetting behavior of a smectic film 
in contact with an attractive wall has been successfully 
studied in Ref. (l^, where the authors found complete or 
partial wetting by smectic depending on the strength of 
the external potential. A infinite (complete wetting) or 
finite (partial wetting) sequence of layering transitions 
was found, some of them ending in a prewetting line. 
Layering transitions at the V-I interface near the V-I-S 
triple point, and thinning transitions in freely-suspended 
smectic films, have been successfully studied using similar 
versions of DFT based on WDA and perturbation theory 

Finally, recent theoretical works have applied related 
models for hard rods in contact with a wall and/or con- 
fined between two walls. These studies were based on 
different approximations : Onsager with restricted orien- 
tations [2l|, [22| , Onsager with Parsons-Lee rescaling and 
free orientations 1231 . |24| and also a WDA functional ap- 
proximation (25l . l26l |. The surface phase diagram of a fluid 
of hard spherocylinders in contact with a single wall pro- 
moting different surface anchoring was analysed in Ref. 
(23| . In Refs. (25l.[26j the surface phase diagram obtained 
for the confined fluid includes capillary nematization and 
smectization of the fluid, and a sequence of layering tran- 
sitions of the confined smectic as the width of the slit pore 
is changed. 

Practically all the experimental work on the wetting 
behavior of liquid crystals has been focused on one- 
component systems, the extension to mixtures being a 
pending issue. Adsorption phenomena in liquid crystal 
mixtures have a fundamental interest since bulk demix- 
ing transitions between two phases, at least one of them 
being smectic, would add much more complexity to the 
surface phase behavior. A recent theoretical work, based 
on Onsager theory, has analysed the phase behavior of 
the I-N interface of binary mixtures of hard spherocylin- 
ders [2Z|. Also, the substrate-isotropic interface of a 
mixture of hard parallelepipeds has been studied within 
the Zwanzig approximation (28j . However, it would be 
interesting to extend these studies to the high-pressure 
regime, where the smectic phase is stable. 

One of the aims of the present work is to elucidate 
the role of the smectic phase in the interfacial phase be- 
havior of binary mixtures. Recent theoretical models of 
mixtures of colloids (spherical or rod-like) and polymers 
[2^43^ (based on the model proposed in Ref. [sj or 
on the recent fundamental-measure functional for hard- 
sphere/hard- needle mixtures jsj) have shown that the 
entropic character of particle interactions, together with 



the coupling between species generated by the external 
surface potential, results in a rich phase behavior. For 
high polymer fugacities, partial wetting of the interface 
between the substrate and the fluid poor in colloidal par- 
ticles by the fluid rich in colloidal particles was obtained. 
In the partial wetting regime, a sequence of up to four 
layering transitions was found. At lower fugacities com- 
plete wetting is reached via a first-order wetting transi- 
tion (located below the critical point). These results were 
confirmed by Monte Carlo (MC) simulations (35| . 

Colloidal rod-like fiuids and their mixtures are paradig- 
matic systems exhibiting liquid- crystal textures similar 
to those of molecular fluids, but the interaction between 
their components have an entropic origin due to short- 
ranged repulsive forces. Intense experimental work on 
pure and mixed suspensions has been done in the last 
two decades, demonstrating this analogy [s^. Also, re- 
cent work has shown the importance of smectic layering 
in the kinetics of the NS phase transition in colloidal 
rods (37I , confirming the analogy between molecular and 
colloidal fiuids as regards the surface-enhanced smectic 
ordering near a bulk phase transition. 

The aim of the present article is to theoretically anal- 
yse the surface phase diagram of a colloidal binary mix- 
ture of rods with the same breath but different lengths 
Li {i = 1,2); in the following we use a length ratio 
s = L2/L1 = 3, with short species being labelled as 1 
and the long species as 2. Particles interact through a 
hard repulsive potential and are constrained to be per- 
fectly aligned along a nematic director, with their main 
axes perpendicular to a hard wall (W), thus simulating 
perfect homeotropic anchoring. This restriction, which 
considerably simplifies the model, is valid as long as one 
is only interested in the surface phase behavior of parti- 
cles exhibiting a high degree of orientational order. The 
study concerns the wetting properties of these mixtures 
when a smectic film partially or completely wets the 
WN interface. Our theoretical tool is based on density- 
functional theory, more specifically on a recently pro- 
posed Fundamental-Measure Functional (FMF) for bi- 
nary mixtures of parallel hard cylinders |38| . 

The impact of restricted orientations was analysed by 
Shundyak and van Roij in the context of the Onsager 
theory [s^ . using the Zwanzig model (discrete orienta- 
tions). It was found to lead to spurious nematic phases 
with very high orientational order. Smectic phases were 
not analysed by Shundyak and van Roij but spurious 
smectic phases might well exist in Onsager theory. The 
parallel particle approximation implicit in the FMF ap- 
proach is not expected to lead to any such anomalous 
phases, since the FMF theory contains a much better 
treatment of correlations, hence of the ordered phases. 

As we will see later, the surface phase diagram of the 
model exhibits three different wetting behaviors depend- 
ing on the value of pressure: (i) At high pressure p we 
find complete wetting by smectic via an infinite sequence 
of layering transitions as the nematic binodal of the bulk 
NS transition is approached. These layering transitions 
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end in corresponding surface critical points characterised 
by values of critical pressure p^'\ n = 1,2, .... (ii) For 

(n) 

sufficiently low pressure such that p < pc , Vn, wetting 
by the smectic film becomes continuous, with adsorption 
coefficients diverging logarithmically. And (iii) for pres- 
sures below the tricritical point, where the bulk NS tran- 
sition changes from first to second order, we find critical 
adsorption by smectic. In this case a modified adsorp- 
tion coefficient diverges logarithmically on approaching 
the second-order bulk NS transition. This divergence is 
a direct consequence of the NS bulk correlation length 
diverging at the transition. 

The article is organized as follows. Sec. [IT] is devoted 
to the presentation of the theoretical model and the nu- 
merical details relevant for the calculation of the bulk 



(Sec. Ill Ap and surface phase diagrams (Sec. Ill Bp . In 
Sec. mil we present the results obtained from numerical 
functional minimization. This section is divided into Sec. 
nil Al where the phase behavior of this particular mixture 
is described, Sec. IIII Bl which contains a description of 
the layering transitions, and Sec. IIII CI devoted to the 
study of the wetting behavior. Finally some conclusions 
are drawn in Sec. IIVI Two appendices are included which 
contain mathematical details on the bifurcation analy- 
sis (Appendix |X| and the derivation of the interfacial 
Gibbs-Duhem relation with composition and pressure as 
independent variables (Appendix [B]) . 



II. THEORETICAL MODEL 



Our particle model consists of a binary mixture of par- 
allel hard cylinders, with both species having the same 
diameter, Di = D2 = D, chosen so as to set the ratio 
of transverse particle area and cylinder length squared 
of the short species to unity, i.e. iiD^ /AL\ = 1. This 
implies a particle aspect ratio of Li/Di = 0.89. Since 
we choose a length ratio s = L2/L1 = 3, the aspect 
ratio of the other particle is L2/-D2 = 2.66. As density- 
functional theory and simulations show, a one-component 
fluid of parallel hard cylinders presents a phase sequence 
nematic-smectic-crystal, which is independent of the as- 
pect ratio. The smectic phase of freely rotating hard 
spherocylinders is known to begin for aspect ratios si 4.1, 
and we should expect a similar behavior for freely rotat- 
ing hard cylinders. Since the phase behavior of a binary 
mixture of parallel particles with identical diameters but 
different lengths depends only on the ratio L2/L1, our 
model might describe a freely rotating binary mixture of 
cylinders with aspect ratios Li/Di > 4.1 and L2/L1 = 3, 
both of which would have a smectic phase. Therefore, 
our choice guarantees that, in the one-component limits, 
the mixture would possess stable smectic phases at high 
enough pressure in the freely-rotating case. 



A. Bulk smectic phase 

A fundamental-measure density-functional theory for 
binary mixtures, in the version proposed in j38| and 
tested against MC simulations in [40|, will be used in all 
calculations. We will consider a mixture which presents a 
non-uniform structure along the z direction. The excess 
free-energy density reads 



<i>Ll = n{- In (1 - 77) + 



Jv_ 

1 - ^ ^ (1 - vf i ' 



(1) 



where we drop the z-dependence for the sake of conve- 
nience and have defined the weighted densities 

= lj2^P*{z^^j2) + p*iz + Kj2)], (2) 

i 

^i') = E r^''\n^')dz', (3) 

J Jz-K,/2 

with T]{z) the local packing fraction of the mixture. In- 
dex i in all sums is assumed to run for i = 1, 2. We have 
defined the dimensionless densities p*{z) = pi{z)L\^ and 
z coordinates are also in units of Li. The Kj parame- 
ter is the particle length of species i in the same units. 
Our choice for Li gives ki = 1 and K2 = s = 3. The 
free-energy functional per unit area can be calculated as 
PJ/A = J dz[$id(z)-f $(z)], with = ksT the in- 
verse thermal energy and 



$id(z) = ^A(2)[lnV,p,(2)-l] 



(4) 



the ideal part of the free-energy density, where Vi is the 
thermal volume of species i. Now we specify for the smec- 
tic phase, which is the lowest symmetry phase considered 
in this work and has the property pi{z + kd) = Pi{z) 
(with d the smectic period and fc G Z). The pressure of 
the mixture can be calculated as 



l3pL\ = d 




^ 3n(z)?7(z) 

ri{z) [l-Ti{z)Y 

2n{z)'q{zf 1 



+ - 



dz. (5) 



During the numerical minimization we have used the fol- 
lowing constraints: (i) the value of the pressure p is fixed, 
and (ii) the composition of the mixture, x = xi ^ Pi/ p, 
is also set in advance. Here p = pi + p2 is the total 
mean density (calculated from the constant-pressure con- 
straint), while Pi ~ d~^ pi{z)dz is the mean density of 
the i-th species. The Gibbs free-energy functional per 
particle, defined as 



P9[Pl,P2] = P" 



[$id(z) + $(2)] rfz + /3pL (6) 
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has been minimized with respect to the densities Pi{z). 
We do this numerically by first discretising the densi- 
ties, defining a grid with points — ~\~ fcA (fc = 
0, ...,A^), and then minimising the function 5(^1,^2) 
with respect to the components of the vectors = 
[pi(zo), • • • ,pi(zAr)], and also with respect to d, using a 
conjugate-gradient algorithm. N is the number of grid 
intervals. The width of the intervals was taken to be 
A/Li ~ 0.01, and A^A — md, where m is the number 
of smectic periods within the minimization box. Varying 
X between and 1 and using the common-tangent con- 
struction for the function /3g{x), we have calculated the 
coexistence values for x and p. Repeating the above pro- 
cedure for different pressures, we obtained the demixing 
binodals. 

When the NS transition is of the second order, one can 
use a bifurcation analysis to find the total packing frac- 
tion 77 = ^ • p*Ki and the smectic period d at bifurcation 
(the local and total packing fractions are equal in the 
nematic phase. In the smectic phase the average of the 
local fraction r](z) over one period gives the total packing 
fraction 77). Also, with the aim to check the relative sta- 
bility of the S with respect to the columnar (C) phase, we 
extended the bifurcation analysis to include the colum- 
nar symmetry. For this purpose we need to solve the 
following set of equations: 



(7) 



where the wave vectors q — (0, q) and q (q_L, 0) are ap- 
propriate for the S and C symmetries, respectively. These 
equations have to be solved for the absolute minimum of 
rf) = det[i/(q, 77)] as a function of q, with i/(q, rj) a 
2x2 matrix defined by the elements 

'1 - picii(q,7?) -pici2(q,77) \ 

i?(q,^)=| , (8) 

-P2ci2(q,?7) 1 - P2C22(q,'7)/ 

with Cij{q^,ri) the Fourier transforms of the direct cor- 
relation functions calculated from the second functional 
derivatives of the free energy functional /3J-[{pi}] with re- 
spect to pi{r) and pj{r'). Expressions for these functions 
and explicit results for the NS and NC spinodals can be 
found in Appendix El 



B. Wall- fluid interface 

The aim is to calculate the equilibrium density profiles 
of the two species in the presence of a hard wall. The 
wall is located at z = and the long axes of cylinders 
are perpendicular to the wall. Thus, perfect homeotropic 
alignment of the nematic director is assumed. This model 
may apply to experimental systems where homeotropic 
anchoring is forced by surface treatment (4l]443 | or by the 
application of an external field (see (45| for an example 
on colloidal discs). 



The values of the chemical potentials of the two compo- 
nents. Pi, will be fixed, which means that the conditions 
of the bulk fluid, far from the wall, will be specified in ad- 
vance and maintained fixed during the minimisation. We 
minimise the grand potential functional per unit area. 



A 



A 



J [vi{z) - pi]pi{z)dz, (9) 



with respect to the density profiles Pi{z). The external 
potentials are defined by 



00, z < Li/2, 
Pv^iz) = { i^l,2. 
0, z> L,/2, 



(10) 



To numerically implement the minimization we proceed 
by first choosing values for the pressure p and the com- 
position of the mixture at bulk, x, and from here calcu- 
lating the chemical potentials pi and the dimensionless 
total density p* = pL\ at an infinite distance from the 
wall, using the following expressions, which apply to the 
bulk nematic phase: 



and 



PPi 



In Xi + In 



P*(l + ^) 
(1 - ,7)3 ' 



(11) 



,1-^ 
p*(4-3?? + ?7^) 
(1 - 



^(3 - 27?) 

(1 - 1? 



(12) 



The implicit Eqn. pT|) has to be solved iteratively to 
obtain p* . In the minimisation the usual boundary con- 
ditions at a large distance H from the wall, Pi{H) ~ Xip, 
have to be imposed. H, the width of the minimiza- 
tion box, is chosen in such a way as to guarantee that 
the structure of the WN interface can be accommodated 
within the box and at the same time the boundary con- 
ditions are satisfied. Finally, the surface tension of the 
interface is calculated as 7wn = 

^[{pl''^}]/A+pH, with 
pI^^ the equilibrium density profiles. 

One of our aims is to obtain the wetting behaviour of 
the mixture when nematic conditions are fixed at bulk 
and the NS demixing transition is approached. This 
means that we need to calculate the surface tension of 
the WS interface for values of the chemical potentials pi 
corresponding to NS coexistence. Therefore pi can be 
calculated from Eqn. (fT2)) . However, if the bulk phase is 
a smectic, and consequently the density profiles are not 
uniform in bulk, the boundary conditions depend on the 
particular value of H chosen, which considerably compli- 
cates the numerical minimization. To avoid this problem, 
we choose to define a symmetric box by using the follow- 
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ing external potentials: 



oo, 



z < L,/2, 



0, L,/2<z<H -L,/2, 
oo, z>H- LJ2. 



(13) 



We minimize fi[{/9i}] with respect to pi{z) by choosing H 
large enough to accommodate two WS interfaces. How- 
ever, due to long-ranged commensuration effects gener- 
ated by the confinement of a layered phase with period 
d in a slit of width H, the minimized grand potential 
exhibits an oscillatory behavior as a function of H, with 
an asymptotically decaying amplitude. To overcome this 
problem, we defined the curve obtained from the local 
minima of and extrapolated to ^ cxd to 

obtain the value of 27^3 (i-e. two times the surface ten- 
sion of the WS interface) . 

To find the surface tension of the NS interface we fol- 
lowed a similar procedure: we defined a box of width H 
with boundary conditions Pi(0) = Pi{H) = p\'^'' (the den- 
sities of the bulk nematic phase coexisting with smectic) 
at both ends of the box. Choosing an initial guess for 
Pi{z), < z < H (close to the profiles of the coexisting 
bulk smectic phase) , we minimized the grand potential to 
obtain 27ns (i-e. two times the surface tension of the NS 
interface). Again the value of H has to be large enough 
to accommodate two NS interfaces. Having the surface 
tensions of all the three different interfaces, one can study 
the wetting behavior of the system, which is discussed in 

Sec. Unci 

Adsorption coefficients will also be used as a convenient 
measure of the wetting and adsorption properties of the 
WN interface. The adsorption coefficients of both species 
are defined as 



r,; 



Pt{z) - P 



(N) 



dz. 



1,2. 



(14) 



In Appendix [B] a derivation is presented of the interfacial 
Gibbs-Duhem relation expressed in terms of the indepen- 
dent variables x and p. Using this equation, a relation 
between the derivative of 7wn with respect to the com- 
position variable x and the adsorption coefficients can be 
obtained: 
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dx 



U{x,p) 



To 



1 



£1 

X 



(15) 



where U{x,p), a function of bulk composition and pres- 
sure, is always positive if the binary mixture is stable 
against NN demixing. This relation has been tested (see 
Appendix |B|) to check for consistency of our numerical 
minimization procedure. Also, the sum rule relating the 
bulk pressure with the densities at the wall (contact the- 
orem), /3p = pi{Li/2) + p2{L2/2), which is automatically 
satisfied by the functional, provides another check for nu- 
merical accuracy. For example, for a mixture with bulk 



pressure fSpLf = 1 and composition x = 0.82, we obtain 
(pi(Li/2) + p2{L2/2))Ll = 0.994, 0.997 and 0.999 for 
values of the discretisation interval along the z axis of 
AL-i = 0.0100, 0.0050 and 0.0025, respectively (obvi- 
ously, in the limit where Az — >■ 0, the sum rule becomes 
exact). 



III. RESULTS 

This section is devoted to the presentation of the re- 
sults obtained from our theoretical model. It is divided 
into three different sections. In Sec. IIII Al we present 
and describe the main features of the bulk phase dia- 
gram. Sees. IIII Bl and IIII CI are devoted to the layering 
transitions and to the wetting behavior, respectively. 



A. Bulk phase diagram 

The bulk phase diagram of the binary mixture, shown 
in Fig. [ll has been calculated using bifurcation analy- 
sis and density-functional minimization, as described in 
Sec. lira Two NS spinodals (dashed curves in Fig. [Ij, 
calculated from the bifurcation analysis, depart from the 
one-component limits 2; = and x = 1 (where, as defined 
above, x is the composition of the mixture as given by 
the fraction of short particles). The values of pressure in 
both spinodals increase as the composition of the mix- 
ture becomes farther from these limits, indicating that 
the two species cannot be easily accommodated into a 
smectic arrangement. The spinodal lines end in a tri- 
critical point (filled circle) and a critical end point (filled 
square), respectively. Functional minimization indicates 
that the Gibbs free energy of the smectic phase is al- 
ways a convex function of composition x in the neigh- 
bourhood of (and above) these lines, which proves that 
the NS transition is of second order, with the smectic 
order parameter increasing from zero at the bifurcation. 
For pressures above the tricritical point but below the 
critical end point, the mixture segregates into a smectic 
phase rich in the long species and a nematic phase rich 
in the short species. 

Two different smectic phases occur in the region of 
smectic stability. These phases are distinguished by the 
relative location of the density peaks of the two species. 
In the smectic phase labelled as Si the profiles are in 
phase, with density peaks of the two species located at 
the same positions, which define the location of the smec- 
tic layers. In the phase called S2 the density profiles 
are out of phase, forming alternatiiig smectic layers: this 
phase exhibits microfractionation |46l - M8j . The smectic 
mixtures with a higher fraction of species 1 (the short 
component) are denoted with a prime in Fig. [TJ Exam- 
ples of density profiles corresponding to these two smec- 
tic phases are shown in Figs. [21[a)-(d). Several regions 
of smectic coexistence exist in the mixture: S1-S2 in a 
narrow pressure interval between a critical and a triple 
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FIG. 1: (Colour online). Bulk phase diagram of the binary 
mixture of parallel hard cylinders in the reduced pressure 
/3pLf-composition x plane (with x — x\, the fraction of short 
particles). Dashed curves represent second-order NS tran- 
sitions, while solid lines are the binodals of the NS or SS 
demixing transitions. The shaded area is the region of in- 
stability. Stability regions of nematic and different smectic 
phases are denoted by letters (see text). Colour lines indicate 
the conjectured wetting behaviour along the NS lines: critical 
adsorption (red), complete wetting by S2 smectic phase with- 
out layering transitions (blue), and complete wetting by Si 
phase mediated by layering transitions (green). Circles: crit- 
ical points. Square: critical end-point. Triangle: tricritical 
point. 



point, S']^-S2 in a corresponding interval between critical 
and triple points, and Si-S'^ and S1-S2 coexistences at 
high pressure (triple points have been indicated by hori- 
zontal dashed lines in the figure). 

We now comment on the possible stability of the 
columnar phase. A complete calculation of the stabil- 
ity of the columnar phase by free-energy minimisation is, 
at present, a highly difficult task. The difficulties stem 
from the computation of two-particle weighted densities 
(38| . Therefore, we have implemented a bifurcation anal- 
ysis, which gives the conditions under which the nematic 
phase becomes unstable with respect to columnar-like 
fluctuations. As shown in Appendix |3 the NC spinodal, 
signalling the instability of the nematic phase against 
columnar-like fluctuations, is always located above the 
NS spinodal for all values of composition. This is an in- 
dication that at least part of the phases depicted in the 
phase diagram of Fig. [T]could be stable, and that the sur- 
face behaviour to be described below could represent the 
real surface behaviour of the model. However, one has to 
be cautious, since a first-order nematic-columnar and/or 
smectic-columnar phase transition could be greatly dis- 
placed with respect to the spinodal lines. Experimen- 
tally, rod-like colloidal particles always have some degree 
of polydispersity. Diameter polydispersity would tend to 
destabilise the columnar phase against the smectic phase, 
while breadth polydispersity would have the opposite ef- 
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FIG. 2: Density profiles of (a) Si, (b) S2, (c) S'a (c), and (d) 
S'l phases in one smectic period. Values of smectic period 
are: (a) d/U = 3.495, (b) 3.835, (c) 1.433, and (d) 1.237. 
In all figures solid and dashed curves correspond to species 
2 and 1, respectively. Values of reduced pressure and com- 
position, {l3pLl,x), are: (a) (1.20,0.06), (b) (1.04,0.28), (c) 
(1.75,0.83), and (d) (1.75,0.88). 



feet [49| . The final balance may depend on several effects 
in a delicate manner. Therefore, one has to be cautious 
until the following aspects are considered: (i) particle 
polydispersity in length and breadth, and (ii) full mini- 
mization of the density functional with respect to density 
profiles having columnar symmetry. We do not pursue 
this analysis here, which is left for future work. 



B. Layering transitions 

In this section, which constitutes the cornerstone of the 
present work, we present a detailed study of the layering 
transitions in the mixture. The stable bulk phase (in 
the region infinitely away from the wall) will be chosen 
to be a nematic phase, characterised by particular val- 
ues of pressure and composition. We first consider the 
case where the pressure is given a value PpL^ = 1.30 and 
the bulk composition x is decreased from a high value 
close to unity. As the nematic branch of the NSi bin- 
odal is approached, a sequence of layering transitions is 
found. At each of these transitions a new smectic layer, 
mostly composed of particles of the long species, appears 
through a first-order (interfacial) phase transition. This 
is illustrated in Fig. [31 where four equilibrium WN in- 
terfaces containing 0, 1, 2 and 3 smectic layers composed 
(essentially) of particles of species 2 are shown. These 
structures will be denoted by WN^ , with i the number of 
adsorbed layers. 

The structure of the WNq interface is interesting. 
Right at the wall there is a mixture of highly localised 
long and short particles with similar densities. For in- 
creasing distance from the wall the density structure be- 
comes much weaker [see Fig. [SJa)]. At the first (WNq- 
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FIG. 3: Density profiles of species 1 (dashed line) and 2 (solid 
line) for /3pL? = 1.30 and (a) x = 0.6800, (b) 0.6400, (c) 
0.61400 and (d) 0.6120 (coexistence value of composition is 
a;coox = 0.6115). The symbol WNi (i = 1,2,3...) denotes the 
interfacial structure containing i adsorbed layers. 




0.016 



FIG. 4: Surface tension of the wall-nematic interface versus 
composition x — a;coex for reduced pressure j3pL\ = 1.3. Sym- 
bols WNi denote the different branches corresponding to wall- 
nematic interfaces containing i adsorbed smectic layers. The 
inset shows the surface tensions for the WNq-WNi surface 
phase transition. 



WNi) layering transition, the wall becomes fully covered 
by long particles and a single very high density peak ap- 
pears [Fig. [3l^b)]. On further decreasing x, the system ex- 
hibits a sequence of phase transitions, WNi-WNi_|_i, each 
involving the addition of a further highly localised peak 
of the long particles. At x = a;cocx = 0.6115 (the com- 
position of the bulk nematic phase coexisting with the 
Si phase at bulk), the wall is completely wet by the Si 
phase, which means that a macroscopically thick smec- 
tic film (consisting of an essentially infinite number of 
smectic layers) is interposed between the wall and the 
nematic phase. We have found up to 12 layering transi- 
tions as a; — ?> Xcoox with x > Xcocx- Access to higher-order 
layering transitions was not possible within the accuracy 
of our numerical procedure. 

In Fig. 0] the behavior of the WN surface tension 
7 — 7wN, as a function of composition, is shown. The 
location of the WN^-i-WN^ layering transition is ob- 
tained from the intersection of the surface tensions cor- 
responding to the two structures. The surface tension 
of the WNq structure, shown in the inset, is somewhat 
peculiar: just before the WNg-WNi layering transition, 
the surface tension exhibits a maximum. This behavior 
can be explained by resorting to Eqn. (jl5p and not- 
ing that the surface tension slope is exactly zero at x* = 

Fi/(Fi+F2) = Pi/(Pi+P2): where = H''^ j" pi{z)dz, 
i.e. when the composition of the mixture at bulk coin- 
cides with its interfacial value, li x > x* , i.e. when the 
interfacial composition is lower than the bulk value, the 
slope of the surface tension is positive, while the opposite 
occurs for x < x*. 

The adsorption coefficients F^, defined in are plot- 
ted in Fig. [5] as a function of x. As can be seen, F2 
abruptly increases at the layering transitions while Fi 
abruptly decreases, i.e. the WN interface exhibits ad- 



sorption of the long species and desorption of the small 
species. This is the natural interfacial path that con- 
nects a nematic phase located far from the wall, and rich 
in short particles, with a smectic film located next to the 
wall, and rich in long species, as a; — > Xcoox- 

Repeating the same calculations, but at a lower value 
of pressure, /SpLf = 1.15, we find that the first two layer- 
ing transitions disappear, while the higher-order transi- 
tions WNi_i-WNj, with i > 2, remain. Fig. [BUb) shows 
that, although the WNq-WNi and WN1-WN2 transi- 
tions are absent, the adsorption coefficients significantly 
increase in the neighborhood of the transition points cor- 
responding to a higher pressure. This behavior is consis- 
tent with the occurrence of critical points for the WNq- 
WNi and WN1-WN2 transitions at critical pressures in 
the interval 1.15 < (3pLl < 1.30. Fig. [D^a) shows a par- 
tial sequence of layering transitions involving up to 12 
surface layers (the maximum number that our numerical 
scheme can deal with). It is reasonable to think that the 
layering transitions will continue up to the bulk transi- 
tion in infinite number (complete wetting scenario). 

A surface phase diagram that includes the first four 
layering transitions is shown in Fig. [T] The following 
trends can be extracted from the figure: (i) all layering 
transition curves approach the nematic binodal as the 
pressure is increased (for a number of layers > 5 the 
curves are too close to the nematic binodal and are not 
visible in the figure); (ii) the critical points, where layer- 
ing transitions terminate, move to lower pressures as the 
number of layers increases for i > 2. It is interesting to 
note that the critical point of the WNq-WNi transition is 
located below that of the WN1-WN2 transition; this fea- 
ture is related to the strong ordering of the WN interface 
just before the WNq-WNi transition. In any case, lay- 
ering transitions terminate at pressures where the bulk 
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FIG. 5: Adsorption coefficients of (a) species 1, (b) species 
2, as a function of composition x — Xcoex (in logarithmic 
scale). Dotted lines indicate the transition points between 
two equilibrium interfaces. The value of the reduced pressure 
is PpLi = 1.3. Labels indicate the number of layers of the 
structures involved in each layering transition. 




FIG. 6: (a) Adsorption coefficient of species 2 as a function 
of X — a;cocx (in logarithmic scale). Equilibrium structures are 
represented by continuous curves, while dashed vertical lines 
indicate layering transitions (metastability branches are not 
meant to be complete), (b) A zoom showing the first two 
layering transitions (in this case adsorption of the first two 
layers does not proceed via surface phase transitions). The 
pressure is fixed at jSpLi = 1.15. Labels indicate the number 
of layers of the structures involved in each layering transition. 



NS demixing transition becomes weak or disappears, i.e. 
when f3pLl ~ 1. 

We note that, depending on the wetting scenario for 
the WS2 interface, the layering transition curves could 
or could not continue above the NS2 spinodal; for exam- 
ple, the wetting regime could change to a partial wetting 
regime, similar to that found in Ref. [l^. Since interfa- 
cial calculations with a bulk smectic phase are difficult, 
we have not carried out this programme in the present 
work. 

Next we briefly discuss the transition strength along 
the layering transition curves. In Fig. [5] the gap in the 
adsorption coefficient of species 2 at coexistence of the 



WNi_i and WNi structures, Ar["'> = T^""'' 
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FIG. 7: Layering transitions (solid curves) between "WNi-i 
and WNi interfacial structures (with i the number of adsorbed 
layers) in the reduced pressure-composition plane. The crit- 
ical points of the transitions are shown with open circles. 
Dashed line is the NS2 second order transition. 




is plotted as a function of composition along the layering 
transition curves and for various indices i. The general 



X 



FIG. 8: Coexistence gap of the adsorption coefficient of the 
second species, AT 2^ , at the "WN^-i-WNi layering transitions 
as a function of composition and for various indices i. Symbols 
refer to layering transitions WNq-WNi (circles), WN1-WN2 
(stars), ■WN2-WN3 (squares) and WN3-'WN4 (triangles). 



trend observed is that, as more layers get involved, the 
transition becomes stronger (i.e. the gap at coexistence is 
larger) . As the index i of the layering transition increases, 
the gap AT 2^ tends to saturate, corresponding to the fact 
that the additional layers adsorbed do not feel the effect 
of the wall and therefore contribute to the adsorption 
coefficient with a constant quantity. 

We end this section with a comment on the origin of the 
layering transitions. As shown in (50| . two-dimensional 
one-component hard-rod fluids in contact with a hard 
wall do not exhibit layering transitions even though the 
bulk transition is of first order [i^. This is also probably 
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the case in the corresponding three-dimensional fluid, al- 
though we have not performed explicit calculations for 
the current model. Therefore, one tentative explanation 
for the phenomenology found in the mixture is that layer- 
ing transitions are due to the entropic coupling between 
the two species mediated by the hard wall: next to the 
wall, where particle densities are high, both species com- 
pete for the gain in entropic volume. Just above the lay- 
ering transition, a mixed layer packs less efRciently, and 
short particles are abruptly depleted from that region, 
with the subsequent abrupt increase in long particles. 
This conclusion would not be affected by the discovery 
of layering transitions in the calculations of Somoza et al. 
[Tgt . who use particles with additional soft, temperature- 
driven interactions; in this fluid the mechanism behind 
the layering transitions would be completely different. 

A look at the structure of the density profiles of the 
WNo interface [Fig. ^'a)] can help us understand this 
effect from a different perspective. The density maxima 
of the two species are clearly separated, due to the dif- 
ferent lengths of the particles (density is maximum ex- 
actly at contact with the wall). However, the maxima 
of the bulk smectic phase are in phase (Si smectic). It 
is only because a sharp change in the interfacial struc- 
ture occurs, via a first-order layering transition, that the 
interface can relax to a structure compatible with that 
in the bulk, i.e. with the correct relative phase. This 
mechanism operates even for structures WN^ with large 
i, when the effect of the wall is not crucial, because the 
density maxima of the two species, in the region between 
the already-formed smectic layers and the nematic, are 
always displaced one with respect to the other. At lower 
pressures, such that the bulk smectic phase is S2, there is 
no such incompatibility between the bulk structure and 
the structure imposed by the wall, and the layering tran- 
sitions vanish. 



C. Wetting behavior 

To obtain a global picture of the wetting behavior of 
the mixture, we have concentrated on four different val- 
ues of reduced pressure: PpLf ~ 1.25 and 1.30 (located 
above the bulk triple point, see Fig. [T|), 1.00 (below the 
triple point and above the tricritical point), and 0.495 
(below the tricritical point) . In the first two cases we have 
found the phenomenology described in Sec. IIII Bl i.e. an 
infinite sequence of layering transitions leading to com- 
plete wetting of the WN interface by the Si phase. Using 
the procedures described in Sec. Ill Bl we have calculated 
the surface tensions of the WN, WSi and NSi interfaces 
for X = Xcocx, which are necessary to discuss the wetting 
behaviour. Their values are collected in Table HI As cor- 
responds to complete wetting by the Si phase, the surface 
tensions fulfill Young's law 7wsi = 7wn + 7nsi (the value 
of 7wN at x = XcocK can be computed by extrapolation of 
7wNi with i ^ 00. In practice i = 32 already gives enough 
accuracy to assess the wetting behaviour. Note that, in 



/3pLf Scocx 


7ws, 


75s, 


7WN 


7wSj + 7NSj 


1.300 0.61150 


0.121166 


0.061370 


0.182535 


0.182536 


1.250 0.58231 


0.122408 


0.055451 


0.177859 


0.177859 


1.000 0.43815 


0.155570 


0.001243 


0.156822 


0.156813 


0.495 0.10000* 


0.181061 




0.181061 





TABLE I: Reduced surface tensions 7* = P'yLf of the WSj, 
WN and NS-, interfaces for different values of the reduced 
pressure. Here j — 1,2 depending on the nature of the smectic 
phase. * indicates value at spinodal. 
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FIG. 9: Density profiles of species 1 (solid line) and 2 (dashed 
line) of the WN interface at bulk pressure /3pLf = 1.00 and 
values of composition (a) Aa; — x — Xcocx ~ 3.2 x lO^'^, (b) 
1.3 X 10"^ (c) 3.2 X lO""* and (d) 6.4 x 10"^ The first and, 
except in (a), the second density peaks are truncated due to 
the small scale of the vertical axis. 

these cases, all the WN^ structures are metastable and 
can be stabilised, even at coexistence, under conditions 
of complete wetting, i.e. when the absolute free-energy 
minimum actually corresponds to i ~ 00). 

The wetting behavior for ppLf = 1.00 is similar to 
that found in (26| and (50l | for one-component hard-rod 
systems: the thickness of the smectic film adsorbed at 
the WN interface grows continuously as a; — > Xcocx and 
diverges at the bulk transition. This behavior is illus- 
trated in Fig. [21 where four density profiles for values of 
composition very close to the bulk transition are shown. 
In this case layering transitions are completely absent; 
adsorption coefficients F^ as a function of x do not ex- 
hibit any discontinuity (Fig. llOp . but diverge logarithmi- 
cally as X — > I'cocx- Young's law for complete wetting is 
also fulfilled within the numerical accuracy that could be 
achieved in this case (see TablelJ). In this case the surface 
tension 7ns 1 is very small and is subject to higher uncer- 
tainties (the value of 7wn was obtained by extrapolation 
to coexistence, x — > Xcocx)- 

Next we discuss the equilibrium density profiles of the 
WN interface at a pressure /3p£f = 0.495 (i.e. below 
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FIG. 10: Adsorption coefficients of species 2 (a) and 1 (b) 
as a function of composition (in logarithmic scale) near the 
wetting transition. The reduced bulk pressure is /^pL'l = 1.0. 



the tricritical point) and, more specifically, the behavior 
of the adsorption coefficients as the bulk NS spinodal is 
approached. Let x* be the composition of the spinodal 
at a given pressure. Since the NS transition is of second 
order, we should find critical adsorption, similar to that 
occurring at the liquid-vapour critical point where the 
adsorption diverges logarithmically as dictated by mean- 
field theory. Far from the wall and close to the bulk spin- 
odal, the WN interface exhibits oscillations with a period 
d* (the smectic period at bifurcation). Thus, the devi- 
ation of density profiles from their bulk values is better 
accounted for by the quantity \pi{z) — pi\ (the analogue 
of the order parameter in the anti-ferromagnetic Ising 
model), and it is convenient to define modified adsorp- 
tion coefficients as 



H 



\Pi{z) - Pi\dz. 



(16) 



The behaviors of r2 and as a function oi x — x* are 
illustrated in Fig. [TT] It can be seen from the figure that, 
while F2 seems to reach a plateau as x — !■ x*, the modi- 
fied coefficient Fg diverges logarithmically as predicted by 
mean-field theory for a critical adsorption phenomenon 
[sH . The value of the plateau is difficult to determine 
due to the huge values of H required to accommodate 
the weakly damped interfacial oscillations that extend 
very far from the wall when x ^ x* . The range of these 
oscillations is of the order of the smectic bulk correla- 
tion length, which diverges ai x = x* . Finally, we have 
checked that the system also exhibits critical adsorption 
in the neighbourhood of the NS' spinodal curve (at higher 
pressure and composition). 

We end this section with a discussion on the impact 
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FIG. 11: The original (a) and modified (b) adsorption coef- 
ficients of second species as a function of bulk composition 
near the NS second order transition at j3pLi = 0.495. In (a) 
the curve is a guide to the eye. In (b) the line is a logarithmic 
fit. 



of the parallel-particle approximation on the wetting be- 
haviour. Complete wetting (either continuous or via a se- 
quence of layering transitions) of a hard wall by a binary 
mixture of hard particles is governed by two factors: (i) 
the effective entropic interactions between particles and 
the wall, and (ii) the distance of the bulk state point from 
the demixing binodal. Since both the free- and restricted- 
orientations models contain these two features, we only 
expect quantitative deviations between the two as far 
as the wetting behaviour is concerned. The situation 
with respect to the critical adsorption phenomenon is 
different, because this is due to the second-order charac- 
ter of the nematic-smectic transition, which may become 
weakly first order for small perturbations of particle ori- 
entations with respect to perfect alignment. In a model 
with free orientations the critical adsorption behaviour 
could be superseded by complete (continuous) wetting. 



IV. CONCLUSIONS 

In this work we have studied the surface adsorption 
phenomena of a liquid-crystal colloidal mixture that has a 
stable smectic phase at moderate pressures. The mixture 
is described by means of a very simple model consisting 
of perfectly aligned hard particles, while the substrate 
is a hard wall inducing perfect homeotropic anchoring 
and nematic boundary conditions far from the wall are 
chosen. Even with these simple assumptions, the bulk 
and surface phase diagrams are so rich that we have con- 
centrated only on the analysis of a single mixture with 
length ratio s = L2/L1 = 3. The theoretical tool used is 
a recently developed fundamental-measure density func- 
tional for mixtures of parallel cylinders (38| . 

We have found a bulk phase diagram with second 
order NS transitions at low pressures, followed by NS 
demixing above a tricritical point. In the low and high 
composition regions of the phase diagram two critical 
points exist, above which two smectic phases, one of 
them micro- fractionated, coexist. Coexistence is ended 
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by corresponding triple points at higher pressures. At 
the highest pressures investigated smectic demixing is 
found, with each smectic rich in one of the species. A 
bifurcation analysis corroborates that the NC spinodal is 
always above the NS spinodal, but does not completely 
clarify the question about the absolute stability of the 
smectic against the columnar phases. In any case, we do 
not expect the NS, and possibly also the SS, demixing 
transitions to be preempted by the columnar phase at 
low pressures. 

The surface phase diagram has three different wetting 
regimes. The first one, located below the tricritical point, 
exhibits critical adsorption as the composition of the bulk 
nematic phase approaches the NS spinodal. In the second 
regime, located approximately above the tricritical point 
and below the triple point (the exact boundaries would 
require further analysis), there exists complete wetting 
of the substrate by a smectic film whose thickness di- 
verges logarithmically as x ~^ Xcocx- Finally, the third 
regime is located above the triple point and is charac- 
terized by the presence of layering transitions that ulti- 
mately lead to complete wetting. A previous theoretical 
study of one-component hard-rod fluids using a density 
functional model [l9| found layering phenomena in the 
semi-infinite system due to strong attractive interactions 
between the wall and the fluid particles. By contrast, 
layering transitions in hard-rod liquid-crystal mixtures 
adsorbed on a hard wall, as shown in the present study, 
is a direct consequence of the wall-mediated entropic in- 
teraction between the two species. 

We expect that the present work serves as a start- 
ing point to initiate experimental studies on the surface 
phase behavior of liquid-crystal colloidal mixtures con- 
sisting of particles that interact through short-ranged re- 
pulsive forces, and having a stable bulk smectic phase. 
These experiments could be guided by the phenomenol- 
ogy found in the present study. Our future theoreti- 
cal studies will analyse the adsorption phenomenology 
of films in the neighbourhood of the bulk triple point, a 
challenging problem that could provide further interest- 
ing phenomena. 
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correlation functions |38l. [53|: 

PiCij{q,v) = {2^ij*o(g±)x {q\\hj) 

+hl,-^i{qi_)x {q\\hl'i) X (9||^j/2)} , (Al) 

where k = Li/ (L), kj = {k + (L) = J^i ^i^i- with 

g|| and q± in units of (L) and D/2, respectively. We have 
defined xi^) = sin(x)/a; and 



Ji(2g) 



+y(l + 2y) 
Ji(2g) 



2j/Jo((7' 



■Ml) 



(A2) 



2(l + 2y)Jo(<z) 



-(l + 6y + 6y2) 



^1(9) 



(A3) 



with y = 77/(1 — 77) and J„(x) the n-th order Bessel func- 
tion of first kind. The NS spinodal can be obtained by 
solving Eqs. ([7]) with = 0, = q, and 

H((Z, 77) = 1 + 512;. {2*o(0);a(9^.) 

i 

+ ^mih\qh/2)} - [AvI/o(0)x(gA//2)]2 , (A4) 

where Al = li — I2, = (P) — 1 being the the poly- 
dispersity coefficient, with {P) = J^i^i^f- ^'^^ columnar 



symmetry g|| =0, \q± 

n{q,v) = 1 + 2*0(9) 



= q, and we find 

*l(g)+A2 [*i(g)-*2(^)] 



(A5) 



To search for a possible NN demixing scenarios, we need 
to solve H(0, 77) = for A as a function of 77, which 
provides the NN spinodal: 



A^ = 



l + 21'o(0) + *i(0) 
*§(0)-4'i(0) 



= (i-l 
77 



1 + 477 + 77^ 
7 — 277 — 77^ 



(A6) 



For our particular mixture {L2 = 3Li), the pressure of 
the NN critical point is PpLf « 10. This result shows 
the metastable character of the NN demixing against NS 
demixing (see Sec. IIII Ap . 

In Fig. [T2] the NS and NC spinodals, in the pressure- 
composition plane, are shown. They never intersect and 
the former is always below the latter. This fact indicates 
that in the low-pressure region of the phase diagram, cal- 
culated in Sec. IIII Al the nematic and smectic phases 
could actually be stable against the columnar phase. 



Appendix A: Bifurcation analysis 



Appendix B: Interfacial Gibbs-Duhem relation 



To implement the bifurcation analysis, we use the fol- 
lowing expressions for the Fourier transforms of the direct 



In this section we obtain the thermodynamic relation 
involving the derivative of the WN surface tension 
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FIG. 12: NS (solid) and NC (dashed) spinodals in the reduced 
pressure-composition plane. 



with respect to composition and the adsorption coeffi- 
cients Ti. Starting from the interfacial Gibbs-Dulrem re- 
lation d'y — — Tid^i, we Irave 



d'j 
dx 



djjLi 
dx 



(Bl) 



For fixed pressure the chemical potentials fii{x,p{x)) 
are only functions of x because p{x) is defined implicitly 
through the constraint p{x, pi^)) ~ Po- Then: 



dp dp/dx 
dx dp/ dp 

Taking into account that 



(B2) 



djii dpi ^ dpi dp 
dx dx dp dx ' 

and using Eqs. ^ and (|B2|), we finally find 

dp, (-l)'-i Alr,^{7^2Tj-,f){l-h) 



dx 



(l-f-47?-t- 772)(l -7y)2 



(B3) 



(B4) 



which, after insertion in (jBip . gives the final result ([T5|. 
with 



U{x,p) = 1 



ff-{7 -2-^-7^) IS? 
(1 + 4r/ + 772)(i _ 



(B5) 



The function rj{x) in Eq. (|B5p can be found from the 
constant-pressure constraint. Comparing Eqs. (jA6p and 
(|B5|) . which contain the same and only density factor 
that can change sign, we conclude that lA{x, P) > only 
when the mixture is stable against NN bulk demixing. To 
check for consistency of our numerical minimization pro- 
cedure, we compare d'y/dx, as calculated from Eq. (jl5p 
[i.e. using the adsorption coefficients Ti obtained from 




FIG. 13: Derivative of the reduced surface tension of the WN 
interface, 7* = PjwnLi, with respect to composition x, ob- 
tained from numerical differentiation of the function 7(1). In 
the inset the difference between this derivative and that ob- 
tained from Eqn. p5[) is plotted (wiggles are due to noise in 
the numerical derivative). The value of the reduced pressure 
is PpLl = 1.3. 



the equilibrium density profiles p[''^] with the numerical 
derivative with respect to x of the surface tension 7 ob- 
tained after minimization. Both results are plotted in 
Fig. [H (for the WN interface) for /SpLf = 1.3. As can 
be seen, both methods reproduce the same function with 
high accuracy, which demonstrates that our calculations 
are fully consistent. Note that the slope of the surface 
tension is equal to zero (and consequently 7 has a max- 
imum as a function of x) for x ~ ri/(ri + r2) [see Eq. 
P^ [. i.e. when the bulk composition is identical to the 
relative fraction of adsorption coefficient of species 1. 
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